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This letter reports results from a search for — > transitions by the MINOS experiment based 
on a 7 x 10 20 protons-on-target exposure. Our observation of 54 candidate v e events in the Far 
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Detector with a background of 49.1 ± 7.0(stat.) ± 2.7(syst.) events predicted by the measurements 
in the Near Detector requires 2 sin 2 (26»i 3 ) sin 2 6» 2 3 < 0.12 (0.20) at the 90% C.L. for the normal 
(inverted) mass hierarchy at Sep = 0. The experiment sets the tightest limits to date on the value 
of 0i3 for nearly all values of Sep for the normal neutrino mass hierarchy and maximal sin 2 (2#23). 

PACS numbers: 14.60.Pq, 14.60.Lm, 29.27.-a 



Observations of neutrinos created in the Sun, in the 
Earth's atmosphere, at nuclear reactors, and by acceler- 
ators provide compelling evidence that neutrinos expe- 
rience quantum mechanical mixing of their weak flavor 
states [IkD]. The resulting neutrino oscillations imply 
that neutrinos have mass and can be represented in ei- 
ther mass or flavor bases, related by the 3x3 PMNS 
neutrino mixing matrix [8|. This matrix parameterizes 
the mixing amplitude using three angles (#12, $23, and 
6*i 3 ), two Majorana phases [§], and a phase (Sep) that 
could give rise to charge-parity (CP) violation in the lep- 
ton sector. The oscillation probability depends on the 
differences in the squared masses of the neutrino states 
and the ratio (L/E) of the distance the neutrino travels 
to the energy of the neutrino. 

In MINOS, the larger mass splitting dominates, and 
oscillations are manifested primarily as the energy de- 
pendent disappearance of muon-neutrinos. MINOS has 
set the most precise measurement of the mass splitting 
|Am 2 | = (2.43±0.13) x l(T 3 eV 2 @, El and requires 
sm 2 (26> 23 ) > 0.9 at the 90% confidence level (C.L.). 
At this mass splitting scale, it is expected that the 
are changing predominantly into v T \ however the sub- 
dominant — > v e transition mode is not excluded [l|. 
Such transitions would indicate a non-zero value of #13, 
the unknown angle of the PMNS matrix and could open 
the possibility of observing CP violation in the lepton 
sector. In this letter, we report new results from the 



search for v u 



v e transitions. 



The most stringent constraint on $13, from the CHOOZ 
reactor experiment 11], implies sin 2 (26*13) < 0.15 at the 
90% C.L. for the value of |Am 2 | measured by MINOS. 
However, a recent global analysis of oscillation measure- 
ments hints at a non-zero value for #13 12j. The CHOOZ 
limit is based on a measurement of the probability for 
electron-antincutrino disappearance. MINOS measures 
the probability of electron-neutrino appearance, which 
additionally depends on sin 2 #23, Sep, and the sign of 
Am 2 . MINOS is the first experiment to probe sin 2 (26*13) 
with sensitivity beyond the CHOOZ limit. An initial 
measurement with 3.14 x 10 20 protons-on-target (POT) 
yielded 35 observed ^ e -likc events with an expected back- 
ground of 27 ± 5(stat.) ± 2(syst.) events 0. This 1.5a 
excess of events is consistent with a value of sin 2 (2#i 3 ) 
near the CHOOZ limit. The present analysis is based on 
an integrated exposure of 7.01 x 10 20 protons-on-target 
and includes the data set from the previous analysis. 

In MINOS, interactions of neutrinos produced in the 
Fermilab NuMI beam line \lM are observed in two detec- 



tors: a Near Detector (ND) with a 29 1 fiducial mass 
1.04 km from the production target and a Far Detec- 
tor (FD) with a 4kt fiducial mass 735 km from the tar- 
get. Both detectors are magnetized tracking calorime- 
ters, composed of planes of 2.54 cm thick steel and 1.0 cm 
thick scintillator (1.4 radiation lengths per plane). The 
scintillator planes are segmented into 4.1cm wide strips 
(1.1 Molicre radii) [l5|. The high statistics data set col- 
lected at the ND establishes the properties of the mostly 
U/j, beam before oscillations. The signature of v ^ — > v e 
oscillations is an excess of v e interactions in the FD rela- 
tive to the expected background based on the ND obser- 
vation. 

Neutrino flavor can be identified in charged current 
(CC) interactions by the event topology produced by the 
associated charged lepton. Muons deposit energy consis- 
tent with a minimum ionizing particle that can be tracked 
through successive detector planes (a track). Electrons, 
on the other hand, deposit energy in a relatively narrow 
and short region (an electromagnetic shower) . Additional 
detector activity can be produced by the breakup of the 
recoil nucleus and other particles produced in the inter- 
action. 

In this analysis, the dominant backgrounds to v e -CC 
events are neutral current (NC) interactions and v^-CC 
interactions with low energy muons. These interactions 
can produce signatures that are similar to those of v e -CC 
events, especially when the hadronic system includes a 
7T°. An irreducible background arises from the 1.3% 
v e +V e component of the beam. This beam v e background 
results primarily from decays of muons produced in pion 
and kaon decays. Their rate below 8 GeV is well con- 
strained by the measured z/ M energy spectrum 0, [lij ]. 
Smaller background components come from cosmogenic 
sources and CC interactions of v T coming from — > v T 
oscillations. 

Selection criteria are applied to events to isolate v e -CC 
interactions and suppress backgrounds. Cosmogenic 
backgrounds in this analysis are reduced to less than 0.3 
events (90% C.L.) in the FD by applying directional re- 
quirements and requiring the events to be in time with 
the accelerator pulse. Selected events must have recon- 
structed energy between 1 and 8 GeV, a reconstructed 
shower, and at least 5 contiguous planes, each with en- 
ergy depositions above half the energy deposited by a 
minimum ionizing particle. Events with long tracks are 
rejected. Further enrichment is achieved using an artifi- 
cial neural network (ANN) with 11 input variables char- 
acterizing the longitudinal and transverse energy depo- 
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sition in the calorimeter [17|. The variables used in the 
ANN are identical to those used in [ill, but the network 
was re-optimized over a sample of simulated events gen- 
erated with a refined detector response model, improved 
event reconstruction, and better modeling of hadron scat- 
tering within the iron nucleus. Maximum sensitivity is 
achieved by selecting events with the neural network out- 
put above 0.7. Background rejection is improved by a 
factor of 1.2 over that reported in [l3| for a similar signal 
efficiency. 

The number of expected background events is deter- 
mined from ND data. The extrapolation to the FD of 
each of the primary background components, v^-CC, NC, 
and beam i/ e -CC, has a different dependence on oscilla- 
tion probability and beam geometry and is treated sep- 
arately Individual background components are deter- 
mined using three beam configurations, each with dif- 
ferent relative background compositions. The first con- 
figuration is the standard one used for the appearance 
search. The hadron production target is located close 
to the first focusing horn, producing a neutrino beam 
peaked at 3 GeV. In the second configuration, the target 
is moved upstream from the horns causing higher energy 
hadrons to be focused and yielding a neutrino spectrum 
peaked at 9 GeV. In a third configuration the current in 
the focusing horns is turned off so no hadrons arc focused. 
Consequently, the low-energy peak of the neutrino energy 
distribution disappears, and the selected event sample is 
dominated by NC events from higher energy neutrino in- 
teractions. 

Data obtained in the above configurations and the sim- 
ulated ratios of rates for each configuration are used to 
extract the three individual background spectra. The 
beam line, detector, and particle propagation simula- 
tion is based on GEANT3 [18| and the hadron production 
yields from the target are based on FLUKA ljjj . Neu- 
trino interactions and further re-interactions of the re- 
sulting hadrons within the nucleus are simulated using 
NEUGEN3 The predicted neutrino energy spectrum 

is adjusted to agree with the ND i^-CC data [Iff. 

Data were collected during three run periods, each with 
somewhat different beam conditions. Most notably, dur- 
ing the third run period, the decay pipe was filled with 
helium at 0.9 atm for safety reasons. The background 
decomposition is performed as a function of neutrino en- 
ergy and is done separately for each run period to account 
for the different beam conditions, a small, gradual target 
degradation, and detector aging. Figure Q] shows the en- 
ergy spectrum measured in the ND for events passing the 
selection criteria and the extracted NC, i^-CC, and beam 
u e -CC components. The ND background is (64±5)% NC, 
(23±5)% ivCC and (13±3)% beam v e -CC events. The 
errors on the components are derived primarily from the 
data and are correlated due to the constraint that the 
background must add up to the observed ND event rate. 
This constraint also leads to a much reduced error on 
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FIG. 1: Reconstructed Near Detector energy spectra of the 
Ve-CC selected events (bold solid). Also shown is the de- 
composition of this spectrum into neutral current (dotted), 
Vfj. charged current (dashed), and beam u e -CC (light solid) 
components determined using the multiple beam configura- 
tion method. The sum of the three background components 
is constrained to agree with the data. Uncertainties on the 
data are statistical and are not visible on this scale; uncer- 
tainties on the components are systematic. 



the FD prediction. A second decomposition technique 
was applied to verify the background components. This 
method uses z/^-CC events with the muon track removed. 
The remnant hits are then processed through the stan- 
dard analysis HI 21 1- This second method yields consis- 
tent ND background components. 

After decomposition of the ND data from each run 
period into separate background components, each spec- 
trum is multiplied by the ratio of FD to ND event rates 
in reconstructed energy bins based on the simulation for 
that component, providing a prediction of the FD spec- 
trum in the absence of v e appearance. Neutrino oscil- 
lations are included when predicting the FD event rate. 
The predictions are summed to give 49.1 expected back- 
ground events, of which 35.8 are NC, 6.3 i^-CC, 5.0 beam 
v e and 2.0 v T (H- 

The efficiency for selecting v e -CC events is estimated 
using remnants from muon-removed i/^-CC events with 
a simulated electron replacing the muon. The embedded 
electron has the same momentum and direction as the re- 
moved muon. Test beam measurements pij ] demonstrate 
that single electrons are well modeled in the MINOS de- 
tectors, and the selection efficiency of electrons agrees 
with the simulation to within 1.6%. With this method, 
we find our efficiency for selecting is e -CC events to be 
(41.6±1.0)%. 

Systematic uncertainties are evaluated by generating 
Monte Carlo (MC) samples in which systematic effects 
are varied over their expected range of uncertainty and 
quantifying the change in the number of predicted back- 
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FIG. 2: Distribution of the ANN selection variable for pre- 
selected events in the Far Detector. Black points show data 
with statistical error bars. The non-shaded histogram shows 
the expected background. The shaded region shows the addi- 
tional v e charged current events allowed from the best fit to 
the oscillation hypothesis as described in the text. 



ground events in the FD. Most of the dominant uncer- 
tainties arise from Far/Near differences. The principal 
systematic effects, listed in Table U include (a) uncer- 
tainties in energy scale, (b) uncertainty in Near to Far 
relative event rate, and (c) uncertainties in the nuclear 
hadronization and intranuclear scattering models. Other 
systematic error sources (d) including uncertainties in 
cross section models, beam flux, and the details of the de- 
tector simulation each contribute to the systematic error 
at lower levels. While uncertainties in the composition 
and kinematic distribution of the particles that emerge 
from the nucleus can be large, these and other uncertain- 
ties associated with neutrino interaction physics mostly 
cancel when comparing the ND and FD data. The use 
of the same materials and detector segmentation in the 
ND and FD is critical in achieving this error cancella- 
tion. The individual systematic errors on the expected 
background are combined in quadrature with the uncer- 
tainty from the decomposition of the background and a 
systematic error on the v T background to give an overall 
systematic uncertainty of 5.6% on the expected number 
of background events in the FD. 

The expected number of background events and its 
uncertainty, along with all the analysis procedures are 
established before examining the full FD data set. Ad- 
ditionally, before counting events in the signal region of 
ANN > 0.7, two FD data samples are examined to check 
the expected event rate and the background rejection in 
the FD. To verify the expected event rate, the full decom- 
position and extrapolation method is applied to events 
well below the signal region (ANN < 0.5), giving a pre- 
diction of 313.6 events. Wc observe 327 events, consistent 
with the prediction to within ltr. Background rejection 
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TABLE I: Systematic uncertainty in the total number of back- 
ground events in the Far Detector. 
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FIG. 3: Reconstructed energy spectrum of the events in the 
Far Detector which pass all v e charged current selection crite- 
ria, with the exception of the energy cut. Black points show 
data with statistical error bars. The non-shaded histogram 
shows the expected background. The shaded region shows 
the additional v e charged current events allowed by the best 
fit to the oscillation hypothesis. 



is verified by examining muon-removed v^-CC events. In 
the FD (92.8 ± 0.9 (stat.))% of muon removed events are 
rejected, in agreement with (93.58 ± 0.05 (stat.))% pre- 
dicted from the ND data. 

Figure[2]shows the number of selected candidate events 
in the FD as a function of the ANN selection vari- 
able. The energy spectrum for the events in the sig- 
nal region is shown in Figure [3l We observe 54 events 
in the signal region with an expected background of 
49.1 ± 7.0(stat.) ± 2.7(syst.), a 0.7er excess over the ex- 
pected background. Similar results were produced by 
a cross check analysis that used a different neural net- 
work based on an alternate event reconstruction algo- 
rithm 24| . Taking into account the improved background 
rejection in the current analysis, this result is consistent 
with the earlier report, based on a smaller data sample. 
From that sample we now select 28 events with an ex- 
pected background of 22.5±4.7(stat.)±l.l(syst.) events. 

Figure @] shows the values of 2 sin 2 (2#i3) sin 2 6*23 and 
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FIG. 4: Values of 2sin 2 (2(9i3) sin 2 #23 and 5cp that produce 
a number of candidate events in the Far Detector consistent 
with the observation for the normal hierarchy (top) and in- 
verted hierarchy (bottom). Black lines show those values 
that best represent our data. Red (blue) regions show the 
90% (68%) C.L. intervals. The CHOOZ limit is drawn for 
Am§ 2 =2.43 x 10 _3 eV 2 , sin 2 (26> 23 )=1.0 

5c p that give a number of events consistent with our ob- 
servation. The oscillation probability is computed using 
a full 3-flavor neutrino mixing framework with matter 
effects (25|, [2^] , which includes a dependence on the neu- 
trino mass hierarchy. Statistical and systematic uncer- 
tainties are included when constructing the confidence 
intervals via the Feldman-Cousins approach (2^. The 
variations of the values of lAm^l, Am^, sin 2 #23, and 
sin 2 (2^12) within their experimental errors [l], [H|-[3] are 
included in the computation of the contours. 

In conclusion, we report improved constraints on 
2 sin 2 (2#i3) sin 2 ^23 from the search for v e appearance by 
the MINOS experiment. The 54 events selected in the 
Far Detector are 0.7 a higher than the expected back- 
ground of 49.1 ± 7.0(stat.) ± 2.7(syst.). Interpreted as an 
upper limit on the probability of — > v e oscillations, 



our data require 2 sin 2 (26>i 3 ) sin 2 6> 23 < 0.12 (0.20) at the 
90% C.L. at 5cp = for the normal (inverted) hierar- 
chy. This measurement represents the best constraint on 
the value of #13 for nearly all values of Sep assuming the 
normal mass hierarchy and maximal sin 2 (2^23) ■ 
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